EXACT FORMULAS FOR TRACES OF SINGULAR 
MODULI OF HIGHER LEVEL MODULAR FUNCTIONS 
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Abstract. Zagier proved that the traces of singular values of the clas- 
sical j-invariant are the Fourier coefficients of a weight | modular form 
and Duke provided a new proof of the result by establishing an exact 
formula for the traces using Niebur's work on a certain class of non- 
holomorphic modular forms. In this short note, by utilizing Niebur's 
work again, we generalize Duke's result to exact formulas for traces of 
singular moduli of higher level modular functions. 



1. Introduction and statement of result 

The classical j-invariant is defined for z in the complex upper half plane 
Hhy 

j( z ) = q- 1 + 744 + 196884(7 + • • • , 

where q = e(z) = e 2mz and J(z) = j(z) — 744 is its normalized Hauptmodul 
for the group T(l) = PSLiffi). All the modular groups discussed in this 
paper are subgroups of T(l). For a positive integer D congruent to or 
3 modulo 4, we denote by Qd the set of positive definite integral binary 
quadratic forms 

Q(x, y) = [a, b, c] = ax 2 + bxy + cy 2 

with discriminant —D = b 2 — Aac. For each Q £ Qd, we let 

-b + is/D 
ZQ= 2a ' 

the corresponding CM point in TC. The group T(l) acts on Qd and we write 
T(1)q for the stabilizer of Q in T(l). The trace of J singular moduli of 
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discriminant — D is denned as 

Qefir>/r(i) 1 wyi 

In [3 Theorem 1], Zagier proved the generating series for the traces of 
singular moduli 

g( z ) : =q- 1 -2- Y t j( D )Q D = ~ 2 + 248 <7 3 - 492 <? 4 + ■ ■ ■ 

D>0 
D=0,3(mod4) 

is a weakly holomorphic modular form of weight 3/2 on To (4). Recently, 
Bruinier, Jenkins, and Ono [2] obtained an explicit formula for the Fourier 
coefficients of g{z) in terms of Kloosterman sums and Duke [3] derived an 
exact formula for tj(D) as follows; 



(1) tj(D) = -24H(D)+ S D (c) sinh 

oo 

c=0(mod4) 



/ 4ttv A D N 



where H(D) = -— ; . — r is the Hurwitz class number and 

v y ^ r i v, 

Sd(c) = ^2 e(2x/c). Using these two results together, Duke [3] 

x 2 =— D(mod c) 

reestablished Zagier 's trace formula Theorem 1]. 

The purpose of this paper is to give a generalization of ([1]) to traces of 
singular values of modular functions of any prime level p. For prime p, 
let Tq(p) be the group generated by Tq(p) and the Fricke involution W p = 
~hp {p ~0 )• Let Qd,p denote the set of quadratic forms Q G Qd such that 
a = (mod p) on which Fg(p) acts. Then the discriminant — D is congruent 
to a square modulo Ap. We choose an integer (3 (mod 2p) with 1 = —D 
(mod 4p) and consider the set Qd,p,p = {[ a >fr> c ] S Qd,p | 6 = /3 (mod 2p)} 
on which Tq(p) acts. For a modular function / for Tq(p), we define the class 
number H p (D) (resp. H*(D)) and the trace tf(D) (resp. tt(D)) by 

Qes D , p /r (p) QeQ D , P /rs(p) 1 ° {P)Ql 
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Here Tq(j))q and Tq(p)q are the stabilizers of Q in To(p) and Tq(p), respec- 
tively. It is easy to see that 



(2) h;{d) 

and 

(3) t* f (D) 



\H P (D), if/3 = 0orp (mod 2p), 
H p (D), otherwise, 



\t f (D), if/3 = 0orp (mod 2p), 
tf(D), otherwise. 



The modularity for traces tf(D) was established by one of the authors in 
0], [5] when Tq(j>) is of genus zero. If / is the Hauptmodul for such Tq(p)* 
and if we define t/(-l) = -1, t/(0) = 2 and t f (D) = for D < -1, then 
the series $^ nr tf(4pn — r 2 )q n ( r , where £ = e(w) for a complex number 
w, is a weak Jacobi form of weight 2 and index p. Meanwhile, using theta 
correspondence, Bruinier and Funke [T] generalized Zagier's trace formula 
to traces of CM values of modular functions of arbitrary level. In particular, 
they showed that if p is an odd prime and f = a{n)q n is a modular 
function for T^(p) with a(0) = 0, then 



£ t}(D) q D + J2(^n)+ P a(n/p))a(-n)- £ 

£>>0 n>l m>ln>l 

-!>=□ (mod4p) 



is a weakly holomorphic modular form of weight 3/2 and level 4p. 

We will obtain in the next section, the following exact formula for t*j(D) 
which is a generalization of JI]). 

Theorem 1.1. Suppose f is a modular function for Tq(jj) with principal 
part Y2m=i a m^(— mz ) a t ioo and define for any positive integers m and c, 

<So(m,c) = e(2mx/c). 

x 2 =— D(modc) 

Then 



f f (D) = £ a. 



N 

"m 



m=l 



c m H*(D) + ^2 S D (m,c)smh 
oo 

c=0(mod4p) 



C 



where c m = —24 ( a(m/p a ) + aim) \ with p a \\m. 

p + 1 
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As an example, consider 

U(37.)J 2 + 37 l V(z) ) ' 

where rj(z) is the Dedekind eta function defined by rj(z) = q2± YinLiO- ~ 1 n )- 
Then / is a modular function for Tq(37) which is of genus 1 and has a Fourier 
expansion of the form q~ 3 — 2q~ 2 — q^ 1 + + 0(q). Since the representatives 
for Qi48,37,o/To(37) are given by [37,0,1] and [74,-74,19], we find from 
equations (J2J) , ([3]) , and Theorem 11.11 that 



38 ^ 

c>0 
c=0(mod!48) 



-SdO-, c) sinh 



sinh 

















2S D (2,c) sinh 



+ / 



8-itVD 



37 + V37f 
74 



where the latter is known to be —2. 

2. Proof of Theorem 11.11 

Throughout this section, T denotes Tq(p). For a positive integer m we 
consider Niebur's Poincare series [6] 

(4) J r m (z,s)= e(-mReMz)(ImMz) 1/2 I s _ 1/2 (2irmImMz), 

where I s -i/2 1S the modified Bessel function of the first kind. Then J- m (z, s) 
converges absolutely for Re s > 1 and satisfies 

(5) T m (Mz, s) = T m (z, s) for M € T and AJ : rn (z, s) = s(l — s)J : m (z, s), 

where A is the hyperbolic Laplacian A = —y 2 (d 2 +d 2 ) for z = x+iy. Niebur 
showed that J- m (z,s) has an analytic continuation to s = 1 [fU Theorem 5] 
and that J- m (z,s) has the following Fourier expansion [6l Theorem 1]; for 
Res > 1, 



(6) T m {z, s) = e(-mx)y 1/2 / s _ 1 / 2 (27rmy) + ^ b n (y, s; -m)e( 



nx 



where b n (y,s; —m) — > (re ^ 0) exponentially as y — » ioo. Hence the pole 
of T m {z, 1) at ioo may occur only in e(— mx)y 1 / 2 / 1 / 2 (2vrmy), which is equal 
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to 

( 7 ) — i/o i/o smh.(2irmy)y 1/2 e(-mx) = . (e(-mz) - e(-mz)). 

So if we multiply J- m (z, 1) by 2-7rm 1//2 , then the coefficient of e(—mz) is nor- 
malized. Now we need to compute the constant term in (2'Km}l 2 )J r m (z, 1). 
It follows from [6, Theorem 1] that bo(y, s, —m) = a m (s)y 1 ~ s / (2s — 1). Here 

(8) a m (s) = 2vr s m s - 1 / 2 m (s)/r( S ) and m (s) = ^ S(m, 0; c)c~ 2s , 

c>0 

where S(m, n; c) is the general Kloostermann sum ^o<d<|c| e (( ma + nd)/c) 
for ("rf) G r. Note that if M = (* G T = rg(p), then M G r (p) or M 
is of the form ( ^{y^ ) with x, y,z,w G Z. In the former case, c is 

v^™ /det=l 

a multiple of p and in the latter case, c = Jpz with p \ z. For n G Z + , let 
u m (n) denote the sum of rre-th powers of primitive n-th roots of unity. We 
observe that 



S(m, 0; c) 



n m (c), if p | c, 

u m (z), if c = A /pz with p\ z. 



J u m (n), if p | re, 

If we define M^n) = < , then 

I p~ s u m (n), if pf re. 



p^ m (s)C(2 S )=p s Y,S(m,0;c)c- 2s £ c'" 
oo c 'ez+ 



2s 



(9) = £ (p s ^(c))c- 2s £ c'- 2s = £ J><(c) ) fc" 

Note that if p f k, then 

fc, if I m, 



2s 



(10) £jA4(<0 = £«*>(<0 = 

c\k c\k 

and if A; = p l k' with Z > 1 and p f fc', then 



0, if k\ m. 



(11) ^fu^c) = ^p s ul l (d) + Y,P S u*M = ^2u m (d) + Y,P S u m (c). 

c\k d\k' c\k d\k' c\k 

p\c p\c 

By adding (p s — 1) Yld\k' u m(d) on both sides of (fTTj) . we obtain 

(p S - l)J2 u m(d) + Y,P Su *m( c ) = £pXn(c). 
d\k' c\k c\k 
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Since u m {d) 
d\k> 

we find that 



k 1 , if k 1 | m, 
0, if k' \ m 



and 



V u m (c) 



c\k 



p s k, if k | m, 
0, if k \ m 



(12) J2p s <^ = < 

c\k 



p s k + (1 — p s )k', if k | m, 
(1 — p s )k' , \i k\m and A;' | m, 

0, if k \ m and k' \ m. 



Writing m = p a m! with p \ m' , we can deduce from (I10p and (I12p that 



A;eZ + c|fe fc'|m' /=! fc'|m' 



2* 



1=1 k'\m' 



2.s 



a^ 2s (m') + (1 -^)a 1 _ 2s (m / )E^" 2S ) / +^ E ^ I ^V^) 1 



-2s 



KKa 



cri_ 2s (m') 
,2s 



1 + (1 - p s 



-2s 



1 - p- 2s 



+ p s ((Ji_ 2s (m) - o-i_2s(m')) 



(13) 



-P 



l+p 



-ai- 2s (m/p a ) +pVi_ 2s (m) 



Recall the constant term in (2irm 1 / 2 )J r m (z, 1) is 

lim 2Trm l / 2 b (y, s, -m) = lim 2-Km 1 / 2 a m (s)y l ~ s / (2s - 1). 

s— >1 S— >1 

By the definition of a m (s) in (JSj) , it is equal to 

lim 2itm 1 l 2 (2-K s m s ~ l / 2 4> m {s)/T(s))y 1 - s /{2s - 1). 

s— »1 

It follows from Q and (fl~3j) that this limit goes to 



47T 2 ? 



m / — p 



■a-i(m/p a ) +pa-i(m] 



K(2) VI 

Thus simple calculations lead us to have the constant term of (2-Km l / 2 )!F m (z, 1) 



as 



(14) 24 
Now we define 



l+p 



■a(m/p a ) + <r(m) 



Cm. • 



•^m^^) = {^m l/2 )T m {z,s) +c m . 
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Then by ©, ([!]), and (J7]), T^iz, 1) is T-invariant harmonic function and 
^-^(z, 1) — e(—mz) has a zero at zoo. Hence it follows from [6l Theorem 6] 
that 

N 



f( z ) = E a m^m( z ^) 



m=l 



for any modular function / for T^(p) with principal part Y2m=l fl m e (" mz ) 
at ioo. Hence 

(15) t}(D) = jr am ( e if?-r^(^.i)]- 

m=l \QeQ D , p /r 1 Ql J 

In order to determine the value — — -^^(zq, 1), we first compute 

for F 
(16) 



QeQ D , P /r 1 
for Re s > 1, 



E ■j=-r = CmH*(D) +27iym V -— ^(zq, a). 

— if) — In 

QeS D , p /r 1 Vl QeQ D , p /r 1 
By Poincare series expansion of f r m (zQ, s) in ([I]), 

27r v / m E -^m^Q, g) _ 27r ^/^ e(-mRe ZQ)(lmzQ) 1 ^ 2 I s _ 1 / 2 (2iTmlmzQ). 
The series in the latter is equal to 

1 /2 

(2mb\(2y/D\ ( 2 v A d\ 

^ e v U 1/2 l 27rm i^ J 

[ap.MeQo.p/rcx, v ^ y V ^ / V ^ / 

= 22 E e 7 -v» 2 ™i^- 

0=1 z(mod2ap) ' ' V ' 7 X 
x' 2 =—d (mod 4ap) 



E -Sbtm.cH 
oo ~ V 



— J 4-1/2 ( 2vrm^^ 



c=0 (mod 4p) 

which converges uniformly for s S [1,2]. Therefore, by (|16|) . 

E ]^n^(^,i) = cm#;p)+ E ^(^.^mh ( ) 

QeQ D , P /r 1 Ql oo V ' J 

c=0 (mod 4p) 



This combined with (|15p completes the proof. 
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